A model of the term structure of interest rates 

based on Levy fields 

Sergio Albeverio 

Institut fur Angewandte Mathematik, Universitat Bonn, Wegelerstr. 6, D-53115 Bonn, 
Germany; SFB 611, Univ. Bonn, Germany; CERFIM (Locarno); Ace. Arch. (USI), 
Switzerland; BiBoS, Univ. Bielefeld, Germany 

E-mail: albeverio@uni-bonn.de 

Eugene Lytvynov 

Institut fur Angewandte Mathematik, Universitat Bonn, Wegelerstr. 6, D-53115 Bonn, 
Germany; SFB 611, Univ. Bonn, Germany; BiBoS, Univ. Bielefeld, Germany 

E-mail: lytvynovOwiener . iam . uni-bonn . de 

Andrea Mahnig 

Institut fur Angewandte Mathematik, Universitat Bonn, Wegelerstr. 6, D-53115 Bonn 
E-mail: andrea.mahnig@gmx.de 



Abstract 

An extension of the Heath-Jarrow-Morton model for the development of instan- 
taneous forward interest rates with deterministic coefficients and Gaussian as 
well as Levy field noise terms is given. In the special case where the Levy field 
is absent, one recovers a model discussed by D.P. Kennedy. 
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1 Introduction 

Heath, Jarrow, and Morton (1992) (see also Heath et al. (1990)) proposed a model of 
interest rates and their associated bond prices in which the price at time s of a bond 
paying one unit at time t > s is given by 



P. 



s.t 



exp 



F s u du 



where F Sjt ,, < s < t, is called the instantaneous forward rate, or just the forward rate. 
Then 

R.:=F a ,„ s > 0, (1.2) 

is called the instantaneous spot rate, or just the spot rate. One also defines the dis- 
counted bond-price process as 



Z S:t :=P Stt exp 



R u du 



o 



(1.3) 



In the Heath-Jarrow-Morton (HJM) model, the forward rates are supposed to 
satisfy the stochastic differential equations 

m 

dF s>t = a(s, t)ds + J2 *) dW l (1-4) 

i=i 

where W 1 , . . . , W m are independent standard Brownian motions and a(s, t) and (3i(s, t) 
are processes adapted to the natural filtration of the Brownian motions. This model 
was, in fact, an extension of the earlier work by Ho et al. (1986). 

Kennedy (1994) (see also Kennedy (1997)), while following the approach of model- 
ing the instantaneous forward rates, considered the case where {F st , 0<s<t<oo} 
is a continuous Gaussian random field which has independent increments in the s- 
direction, that is, in the direction of evolution of 'real' time. This framework includes 
the HJM model in the case where the coefficients a(s,t) and (3i(s,t) in (1.4) are de- 
terministic. An important example of application of the Kennedy model is the case 
where the forward rates are given by F s t = n s t + X s>t with /j, s t being deterministic 
and X Sjt a Brownian sheet (see, e.g., Adler (1981) for this concept). In the latter case, 
F Syt has independent increments also in the t direction. Furthermore, this may be in- 
tuitively thought of as the situation of (1.4) driven by an uncountably infinite number 
of Brownian motion. Kennedy (1994) gave a simple characterization of the discounted 
bond-price process to be a martingale. In particular, he showed that the latter is true 
if and only if the expectation fi Syt of F Syt , < s < t < oo, satisfies a simple relation. 

In Bjork et al. (1997a, b) (see also Bjork et al. (1999)), it was pointed out that, in 
many cases observed empirically, the interest rate trajectories do not look like diffusion 
processes, but rather as diffusions and jumps, or even like pure jump processes. There- 
fore, one needs to introduce a jump part in the description of interest rates. The authors 
of these papers considered the case where the forward rate process {F s t , < s < t} 
is driven by a general marked point process as well as by a Wiener process (Bjork et 
al. (1997b)), or by a rather general Levy process (Bjork et al. (1997a)), and the matu- 
rity time t > is a continuous parameter of the model. In particular, an equivalence 
condition was given for a given probability measure to be a local martingale measure, 
i.e., for the discounted bond-price process {Z s t , < s < t} to be a local martingale 
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for each t > (Propositions 5.3, 5.5 in Bjork et al. (1997a), see also Theorem 3.13 
in Heath et al. (1992)). This condition, formulated in terms of the coefficients for the 
forward rate dynamics, generalizes the result of Heath et al. (1990) which was obtained 
for the diffusion case. 

Other generalizations of the HJM model in which the forward rate process satis- 
fies stochastic differential equations with an infinite number of independent standard 
Brownian motions (i.e., m — oo in (1.4)) were proposed in Yalovenko (1998), Kusuoka 
(2000), Popovici (2001), and Liitkebohmert (2002). 

In the present paper, we follow the approach of Kennedy (1994, 1997), but suppose 
that the forward rates {F Syt , < s < t < oo} are given by a Levy field without 
a diffusion part. In particular, {F Syt } has independent increments in both s and t 
directions. Analogously to Kennedy (1994), we give, in this characterization of 

the martingale measure. We also show that, under a slight additional condition on the 
Levy measure of the field, it is possible to choose the initial term structure {^o,t> t > 0} 
in such a way that the forward interest rates are a.s. non-negative. This, of course, was 
impossible to reach in the framework of the Gaussian model, which caused problems 
in some situations (see Sect. 1 of Kennedy (1997)). We then present two examples 
of application of our results: the cases where F s>t is a "Poisson sheet" (this case was 
discussed in Mahnig (2002)), respectively a "gamma sheet." Finally, we mention the 
possibility of unification of the approaches of Kennedy and of the present paper, by 
considering F Sjt as a sum of a Gaussian field and an independent Levy field, and thus 
having a process with a diffusion part as well as a jump part. 



2 The model based on Levy fields 



Let £>:=Cq°(IR 2 ) denote the space of all real-valued infinitely differentiable functions 
on IR 2 with compact support. We equip T> with the standard nuclear space topology, 
see, e.g., Berezansky et al. (1996). Then T> is densely and continuously embedded into 
the real space L 2 (R 2 , dx dy) . Let V denote the dual space of V with respect to the 
"reference" space L 2 (R 2 , dxdy), i.e., the dual pairing between elements of V and V is 
generated by the scalar product in L 2 (R 2 , dxdy). Thus, we get the standard (Gel'fand) 
triple 

V D L 2 (R 2 , dxdy) D V. 

We denote by (■, ■) the dual pairing between elements of V and V. Let C(V) denote 
the cylinder cr-algebra on T>' . 

We define a centered Levy noise measure as a probability measure v on (V ,C(V)) 
whose Fourier transform is given by 



/ 



v{duo) = exp 



(e iTv( - x ' y) -l-ir(f(x,y))a(dr)dxdy 



^ev (2.1) 
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(see, e.g., Ch. Ill, Sec. 4 in Gel'fand and Vilenkin (1964)). Here, R+ = (0, +oo) and a 
is a positive measure on (R + , £>(R+)), which is usually called the Levy measure of the 
process. We suppose that a satisfies the following condition: 

/ T 2 a(dr) < oo. (2.2) 
Jr + 

The existence of the measure v follows from the Minlos theorem. 
For any ip EV, we easily have 

/ (id,ip) 2 v(dtd) = / t 2 a(dr) / ip(x,y) 2 dx dy. (2.3) 
Jv> Jr+ Jr 2 

Thus, the mapping / : L 2 (R 2 ,dxdy) -> L 2 (V',du), Dom(J) = £>, defined by 

(Icp) (u) :={u>, <f), ip eV, u eV, 

may be extended by continuity to the whole L 2 (R 2 , dx dy). For each / G L 2 (R 2 , dx dy), 
we set (•,/):=//. Thus, the random variable (r.v.) (a;,/) is well-defined for u-a.e. 
id ET>' and equality (2.3) holds for / replacing <p. 

Let k : [0, oo) 2 — > [0, oo) be a measurable function which is locally bounded. For 
each s,t > 0, we define the r.v. X Sjt as follows: 

X S;t (o;):=(a;(a;,y),l[o iS ](x)l[o,t](y)x(a;,y)), i/-a.e. u e V , (2.4) 

where x, y denote the variables in which the dualization is carried out. It follows 
from (2.1) that X s>t is centered and has independent increments in both the s and t 
directions. 

We note that, in the case where x(x, y) = 1, {X Sft , < s < t} is a Levy process 
for each fixed t > 0. Indeed, it follows from (2.1) that the Fourier transform of X s t is 
given by 



[ e iXXs ^ v(did) = exp st [ (e tTX - 1 - ir\) a{dr) 



A e 



In particular, the Levy measure of the process {X s t , < s < t} is equal to £<r. 
Let F Sit be the forward rate for date t at time s, < s < t. We suppose that 

i^,* = fis,t + Xs,t, 0<s<t, (2.5) 

where /i Stt is deterministic and continuous in (s,t) on {(s,t) G R 2 : s < t}. The price 
at time s of a bond paying one unit at time t > s is then given by (1.1). We note that 
the random variable J^' X SiU du is z/-a.s. well-defined and 



I 



X s , a (td)du= (id(x,y),l [0jS] (x)l [0i t}(y)x(x,y)(t-(sVy))) for u-a.e. id E V . (2.6) 
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Indeed, for each / G L 2 (R 2 , dx dy) , we have by (2.3): 
J Qf X a>u (w) du^ (cj, /) u(du) 

X SjU (uj){uj, f) u(dw) du 
= / r 2 a(dr)- / / l [0jS] (x)l [oM (y)x( y x,y)f(x,y)dxdydu 

Jr + Js Jr 2 

= r 2 a(dr)- l[o,s](x) ( / l[o,u](y) ^ ) v) dx dy 

Jr+ Jr 2 \Js J 

= T 2 a(dr)- l [0iS] (x)l m (y)(t - (s V y))x(x,y)f(x,y) dxdy, (2.7) 
Jr + Jr 2 

which implies (2.6). 

Analogously to (2.6), (2.7), we have: 

\ R u du = fJ, u ,udu+ {u(x,y),l[o t3 ](x)l[o >a ](y)x(x,y)(s - (x Vj/)), 
Jo Jo 

where the spot rate R u is defined by (1.2). Thus, the discounted bond-price process 
Z Sit is well-defined by (1.3). 

We denote by T s , s > 0, the a-algebra generated by the r.v.'s F UjV , < u < s, 
u < v, which describes the information available at time s. 

Theorem 2.1 The following statements are equivalent: 

(a) For each t > 0, the discounted bond-price process {Z st , T s , < s < t} is a 

martingale. 

(b) 

Hs,t = lM>,t+l I f Tx(x,y)(l-e~™ (x ^ t - (xV y»)a(dr)dxdy 
Jo Jo J«. + 

for all s, t > 0, s < t. 
(c) 

P Stt = E(e~ £ Ru du | Fs) for all s, t > 0, s < t. 

Proof. We first show the equivalence of (a) and (b). Absolutely analogously to the 
proof of Theorem 1.1 in Kennedy (1994), we conclude that (a) is equivalent to the 
following condition to hold: 

/ exp ( - [ (F suu - F S2yU ) du - [ " F s2,u) du] du = 1 (2.8) 

JV V Js 1 Js 2 J 
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for all < s 2 < si < t. By (2.4) and (2.5), (2.8) is equivalent to 



exp - / (uj(x,y),l [s2jSl] (x)l [0M (y)x(x,y))du 



si 



(w(x,y),l[ aatU ](x)l[ 0tU ](y)x(x,y)) duj u(dto) 
exp ^ jf (ii suu - Hs2,u) du + (u U)U - u S2jU ) duj (2.9) 



for all < s 2 < si < t. Analogously to (2.6) and (2.7), we have 



(u(x, y),l [s2jSl] (x)l [0M (y)x(x, y)) du + (u(x, y),l [s2iU] (x)l [oM (y)x(x, y)) du 

J S2 

= (uj(x,y),l [s2:Sl] (x)l [0jt] (y)x(x,y)(t- (x\/y))) for u-a.e. to E V. 
Hence, it follows from (2.1) that condition (2.9) is equivalent to 



/ / [ \e~ T * {x ' y){t ~ {xS/v)) -l + Tx(x,y)(t-(xVy)))dxdya(dT) 
Jr + Jo Js 2 

= / (V si ,u ~ fJ>s 2 ,u) du + / (a u , a - u S2iU ) du (2.10) 

J 51 J S2 

for all < s 2 < sx < t. We remark that the function under the sign of integral on the 
left hand side of (2.10) is integrable. Indeed, let us set C t := sup^ )2/6 [o,t] x(x,y)- Then, 
for all r E (0, 1] and x, y E [0, t], we have 

\ e -rx( X ,y)(t-(xVy)) _ l + T „( Xj y)( t _ (xVy))\ 

< < r 2 C?t 2 exp(rC t t) < r 2 C 2 t t 2 exp(C t t). (2.11) 



n=2 



This, together with the fact that J, Q 1 , r 2 a{dr) < oo, yields the integrability of the 
integrand on the left hand side of (2.10) on (0,1] x [0,t] x [s!,s 2 ]. Furthermore, for 
r E (1, +oo) and x,y E [0, £], we have 



-Tx{x,y)(t-{x\Jy)) 



l + rx(x,y)(t- (xVy))\ < l + rtC t . 



This, together with the fact that +oo ^ r a(dr) < oo, completes the proof of the 
integrability of the integrand on the left hand side of (2.10). 

We now fix t > and suppose, for a moment, that /i Sjt has the following form: 

Hs,t = Ho,t+ / ^ t (x,y)dxdy, (2.12) 
Jo Jo 
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where ^/ t (x,y) is an integrable function on [0,t] 2 . Then, 

ft psi rt psi 

U,U fJ*S2,U 

)du= / %(x,y)(t-(xVy))dxdy. (2.13) 

</«l J S2 J J S2 

Comparing (2.13) with (2.10), we see that condition (2.10) is, at least formally, satisfied 
if ^t(x, y) has the form 

%(x,y) = (t-(xVy))- 1 [ (e^ x ' v ^ xV ^-l + Tx(x,y)(t-(xVy)))(r(dT). (2.14) 

To show that this inserted into (2.12) gives indeed a solution of (2.10), we have to 
verify that the ty t (x,y) given by (2.14) is integrable on [0,t] 2 . Analogously to (2.11), 
we get 
't rt 



m\(t-(xV y yy\ e -r^y)(t-(xvy)) _ 1 + rx ^ 5 y )(t-(xVy)))\ a(dr) dx dy 

^T n x{x,y) n {t-{xV y)) 71 - 1 . , 

> v yj y , — aidr) dx dy 

>o Jo Jmil n\ 

< t 3 C 2 e tCt [ r 2 a(dr) < oo. (2.15) 
-'foal 




Next 



m\(t- (X \jy)yl( e -™(x,y)(t-( x vy)) -i + T x( Xj y)(t-( x \/y)))\a(dT)dxdy 
a.+°°) 

< I ff \{t- (xV y))-\e-™ {x ' y){t - {xyy)) - l)\a(dr)dxdy 

JO JO i(l,+oo) 

+ t 2 C t [ ra(dr) 



'(l,+oo) 

a, 



< 2t 2 C t / ra(dr), (2.16) 

J(l,+oo) 

where we used the estimate: 1 — e~ a < a for all a > 0. Thus, by (2.12) and (2.14)- 
(2.16), statement (a) holds for 

H.,t = to,t+f f [ (t-ixVyyy^e-™^-^ 
Jo Jo Jr+ 

(2.17) 

Let us now suppose that (a), or equivalently (2.10), holds. We set s 2 = and 
Si = s. Then, (2.10) takes the following form: 

m( e -r*(* lW )(t-(xv«)) _ x + rx (a; >y )(t _ ( x V y))) <r(dr) cfedy 



= / (a* s ,u - A*o,u) du + / (// U)U - Ho, u ) du. (2.18) 
Differentiating (2.18) in t yields for s < t: 

fis,t-to,t= f f [ T^y^l-e-^^-^y^aid^dxdy. (2.19) 
Jo Jo Jm. + 

That the integral on the right hand side of (2.19) is finite may be shown analogously 
to (2.15), (2.16). 

Since for the /i Sji given by (2.17) statement (a) holds, this /i Sjt also satisfies (2.19). 
Therefore, 

/ / / (t- (x Vy))-\e- T ^ x > y)it -( xyy)) - 1 + rx{x,y)(t- (x V y))) a(dr) dx dy 
Jo Jo Jr+ 

= [ [ [ Tx(x,y)(l-e~ T ^ x ' y)it - {xVy)) )a(dT)dxdy, 
Jo Jo Jr + 

which implies the equivalence of (a) and (b). 

Let us now show the equivalence of (b) and (c). Analogously to Kennedy (1994), 
we conclude that (c) is equivalent to 

exp ^ - j (F U;U - F s ^ du = 1 (2.20) 

for all s, t > 0, s < t. Analogously to the above, we see that (2.20) is, in turn, equivalent 
to 



/ / / (e- T ^ x ' y) ^- {xyy)) -l + Tx(x,y)(t-(xVy)))a(dT)dxdy 

Jo J s Jr + 

= J (a*u,« - AO du ( 2 - 21 ) 

for all s,t > 0, s < t. Setting in (2.10) S2 = s and S\ = t, we see that (2.21) is a 
special case of (2.10), so that (b) implies (c). To show the inverse conclusion, we follow 
Kennedy (1994). Differentiating (2.21) in t yields 

f f [ rx(x^)(l-e-^>«)(*-^))) < 7(dr) < tedy = /i t)t -^ (2.22) 
Jo Js Jr + 

for all s, t > 0, s < t. Setting s = in the latter equation gives 

f f [ Tx(x,y)(l-e- T ^ x ^ t -^)a(dT)dxdy = ^ t -^ (2.23) 
Jo Jo Jr + 

Subtracting (2.22) from (2.23) implies (b). □ 



Corollary 2.1 Suppose that the Levy measure a additionally satisfies 

(r) CT := / ra(dr) < oo. (2.24) 
Jr+ 

Suppose that statement (a) of Theorem 2.1 holds and suppose that the initial term 
structure {/xo,t, t > 0} satisfies 



U ,t> [ [ x{x,y)dxdy(r) a , t>0. (2.25) 
Jo Jo 



Then the forward rate process {F st , 0<s<t<oo} and the spot rate process {R t , t > 
0} take on non-negative values v-a.s. 

Proof. By (2.5) and Theorem 2.1, we get 

F 8 , t = A*o,t- f f I rx(x, y) e -™(x,y)(t-(*vy)) a{dr) dx dy+ x %t , t>0,0<s<t, 
Jo Jo Jwi+ 



where 



X s j:=X s j+ / x(x,y)dxdy ■ (r) a . 
Jo Jo 



Under condition (2.24), the measure v is concentrated on the set of all signed 
measures of the form J2^=i T n^(x n ,y n ){dx dy) — (r) a dxdy, where S a denotes the Dirac 
measure with mass at a, r n e supper, n e N, and {(x n , y n )}^Li is a locally finite set in 
M 2 , see, e.g., Lytvynov (2003). Therefore, by (2.4), X s>t takes on non-negative values 
z/-a.s. Furthermore, it follows from (2.25) that 

Ho,t ~ [ [ [ r *( x > y)e- T ^ v){t - {xyy)) a(dr) dxdy>0, t > 0, < s < t, 
Jo Jo Jr+ 

from where the statement follows. □ 

Let us consider two examples of a measure v satisfying the assumptions of Theo- 
rem 2.1 and Corollary 2.1. 

Example 1. (Poisson sheet) We take as v the centered Poisson measure tx z with 
intensity parameter z > 0, see, e.g., Hida (1970). The Levy measure a has now the 
form z8\. Thus, the Fourier transform of n z is given by 



/ e l{ ^) ■K z (duS) = exp / {e iLp{x ' y) - 1 - i(p(x, y))zdx dy 

JV L J M 2 



up e v. 



We set x(x, y) = 1. Then, X s>t given by (2.4) with the underlying probability measure 
v = tc z is, by definition, a Poisson sheet, and for each fixed t > 0, {X Sjt , < s < t} is 
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a centered Poisson process with intensity parameter tz. Statement (b) of Theorem 2.1 
now reads as follows: 

t*8,t = A*o,t + z((2- s)e s ~ t - 2e~ l - s + st). 
Condition (2.25) now means //o,t 

> zt 2 , t > 0. 

Example 2. (Gamma sheet) We take as v the centered gamma measure ^ z with 
intensity parameter z > 0, see, e.g., Lytvynov (2003). The Levy measure a on R + has 
the form 

a(dr) = z dr. 

T 

The Fourier transform of 7 2 may be written as follows: 

/ e l{ul,lp) lz(du) = exp I - / (log(l-^(x, y))+<p(x, y)) z dxdy ) , y? e P, M < 1. 
Jo' V </r 2 / 

We set x(x, y) = 1. Then, X Syt given by (2.4) with the underlying probability measure 
v = 7 2 is, by definition, a gamma sheet, and for each t > 0, < s < £} is a 

centered gamma process with intensity parameter tz. Statement (b) of Theorem 2.1 
now reads as follows: 

l*s,t = fMj,t + z (ri + 2s + 2(1 + t) log ^H" 5 ^ - s log(l + t - s)^j . 

Condition (2.25) means /io^ > zt 2 , t > 0. 

It is possible to construct a model of forward interest rates which unifies the ap- 
proach of Kennedy (1994) to modeling the forward interest rate with our approach. 
Indeed, consider F s>t in the form 

F a ,t = Hs,t + X 8 , t + Y att , 0<s<t, (2.26) 

where /i s>t and X s>t are as in formula (2.5) (thus, as in our approach) and Y s>t is a 
centered continuous Gaussian random field that is independent of X u>v , < u < v < oo, 
and has covariance 

Cov(Y Sl>tl , Y s , 2jt2 ) = c(si A S2, ti,t 2 ), 0<Si<U, i = l,2, 

with a function c satisfying c(0,ti,t2) = (as in Kennedy's approach). 

The following theorem may be easily proved by combining the proof of Theorem 1.1 
in Kennedy (1994) and the proof of Theorem 2.1. 

Theorem 2.2 Theorem 2.1 remains valid for the forward rates {F s t , < s < t < oo} 
given by (2.26) if we set the deterministic term fj, Stt in statement (b) to be 

Hs,t = fM>,t+ III rx(x,y)(l-e^ x ^ t - {x ^)a(dr)dxdy+ I c(sAu,u,t)du 
Jo Jo Jr+ Jo 

for all s, t > 0, s < t. 
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